Experimental NMR implementation of a robust quantum search algorithm 
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Graver's quantum search algorithm, involving a large number of qubits, is highly sensitive to errors in 
the physical implementation of the unitary operators. This poses an intrinsic limitation to the size of the 
database that can be practically searched. The lack of robustness of Graver's algorithm for a large number 



*£h of qubits is due to quite stringent 'phase-matching' condition. To overcome this limitation, Tulsi suggested 



a modified search algorithm [PRA 78, 022332] which succeeds as long as the errors are reproducible and 
reversible while Graver's algorithm fails. Such systematic errors arise often from imperfections in apparatus 
setup e.g. the errors arising from imperfect pulse calibration and offset effect in NMR systems. In this paper, 
we report the experimental NMR implementation of the modified search algorithm and its comparison with 
the original Graver's algorithm. We experimentally validate the theoretical predictions made by Tulsi. 
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I. INTRODUCTION 



Quantum computation has developed as an exciting field of research in the last decade and 
it has generated wide interest among scientists and engineers. It offers the opportunity of cre- 
ation of algorithms that are radically different and more efficient as compared to their classical 
counterparts. Shor's prime factorization algorithm Jit and Graver's quantum search algorithm [2] 
have theoretically demonstrated the power of quantum algorithms. However, the experimental im- 
plementation of the quantum algorithms is still quite challenging. Nuclear Magnetic Resonance 
(NMR) has been the vanguard among the presently available techniques for physical implemen- 
tation of quantum algorithms J4, 5]. Till date, the algorithms have been tested on systems with a 
small number of qubits with a presumption that once a quantum computer with large number of 
qubits are made, more real world application of the algorithms can be implemented. Implementa- 
tion of the quantum algorithms on very large system requires the application of a large number of 
unitary operators. As any physical implementation involves some amount of error which accrue 
when the unitary operators are applied in tandem, physical implementation of an algorithm in a 
large system becomes difficult. The sensitivity of the algorithm to small errors can lead to it's 
failure. 

Graver's quantum search algorithm, or more generally the quantum amplitude amplification 
algorithm, is designed to search a marked item from an unsorted database [Q, 3I. It drives a 
quantum computer from a prepared initial state I*) to a desired target state \t), which encodes the 
marked item. Generally, \s) is prepared by applying a unitary operator U on a particular basis state 
|0), i.e. \s) = U\0). The heart of the algorithm is the Graver's iteration operator Q given by 

g = i s i t = ui u t i t , i^ = i-2\ifr)(ifr\. (i) 

Thus 1^ is the selective phase inversion of the state If a = \(t\s)\ then n/Aa times iteration 
of Q on \s) yield the target state \t) with a high probability. For searching a database of N items, 
the initial state \s) is chosen to be the equal superposition of all basis states each of which has a 
probability amplitude 1 / V^V- It is generated by applying the Walsh-Hadamard transform W on 
the basis state |0), i.e. \s) = W\0). Since \t) is a unique basis state, a = 1/ "vW and 0{ "vW) times 
iterations of Q = WI WI t on \s) yield the target state \t). 

In this paper, we consider the case when the implementation errors cause the deviations in 
selective phase inversions, I s and In other words, we want the apparatus to implement {I s ,I t } 



but due to errors, the apparatus implements {/; , If] where 



i; = l-(l- e m )W)m (2) 

is the selective phase rotation of by angle oj. Then the Graver's operator becomes Q = if If 
and the well-known phase-matching condition [6] demands 

<P - if = 0(a) (3) 

for Graver's algorithm to succeed. For large database size, N » 1 and a = 1/ "vW «: 1 and the 
above condition becomes quite stringent. From the implementation point of view, satisfying Eq. [3] 
is tough as the phase rotations on state \s) and \t) are not equal in general. Therefore, as the size 
of the database increases, there is a high risk that Graver's algorithm fails even if there are very 
small errors in the implementation of the operators. 

To take into account the above mentioned problem, Tulsi has modified the quantum search 
algorithm The algorithm is based on the assumption that errors are (i) reproducible and (ii) 
reversible. The reproducibility allows us to implement the transformations [if, If} at our disposal 
while the reversibility allows us to implement the inverse transformations {77 , If v ) at our disposal. 
Then the collective effect of the errors can be cancelled by iterating the following operator 

t = r/i^ifTf = ui~ 1/ V ui* tfif . (4) 

Note that for (f> = <p = n, T is just two steps of Graver's algorithm, i.e. T = Q 1 . Tulsi has shown 
that n/4a sin | sin | times iteration of T on |s) yield the target state \t) with high probability. 
Therefore, if a is small (i.e. the database is large), small difference between (p and (p can cause 
the Graver's algorithm to fail while the modified algorithm still succeeds in finding the target state 



(see Fig. |l(a)| for simulation results). However it may be pointed out that Graver's algorithm is 
self correcting if (f> = (p (Fig. |l(b)[ ). The complexity of both the algorithms remains almost the 
same for {<f>, <p) n QV|] . 

It should be noted that for the experimental demonstration of the difference between the orig- 
inal and the modified search algorithm for large database, it is not necessary to implement them 
on a very large system. We can simulate the behaviour of the algorithms for large database by 
preparing a small system with a = \(t\s)\ «; 1, i.e. initially, the target state has a low probability 
amplitude. As there are no other restrictions on the size of system, a two qubit system is suitable 
enough for this purpose. 
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FIG. 1: Simulation results for the original and the modified search algorithm. The line joining the points is 
just a guide for the eye. (a) Here a = 0.00091, <p = n and <p = 0.9k. Thus the phase-matching condition 
<p - tp - O(a) is not met. For the modified algorithm, the probability of the marked state reaches very near 
to 1 whereas for the original algorithm, it is always low. (b) Here a = 0.00091 and (f> - (p - 0.9n. Thus the 
phase-matching condition is satisfied. Graver's algorithm is self correcting even for large errors as long as 
the phase matching condition is met and its performance is quite similar to that of the modified algorithm. 

To experimentally verify the algorithm of Tulsi fl7J], the original and the modified search 
algorithms are implemented here in an NMR Quantum Information Processor. The implementa- 
tion procedure consists of (i) preparation of the pseudo-pure state (PPS), (ii) preparation of the 
superposition of all the states such that the marked state has a low probability amplitude, (iii) 
application of the original/modified iterations and finally (iv) measurement. The experiment has 
been carried out at room temperature in 11.7 Tesla field in a Bruker AV500 spectrometer using a 
QXI probe. The system chosen for the implementation of the algorithm is Carbon- 13 labeled chlo- 
roform ( 13 CHC1 3 ), where and 13 C act as the two qubits. The ! H and 13 C resonance frequencies 
at this field are 500 MHz and 125 MHz respectively and the scalar coupling between the spins is 
Jhc= 209 Hz. The NMR Hamiltonian for a 2-qubit weakly coupled spin system is [8], 



M 3 = vj] + v 2 I 2 7 + J n l\%, 



(5) 



where v, are the Larmour frequencies and the J 12 is the scalar coupling. The equillibrium density 
matrix, which is the starting point of any algorithm in NMR quantum information processor, under 
high temperature and high field approximation is in a highly mixed state represented by [8D, 



Peg = JhI H z + JcI C z = 7h (/f + 0.25/f ) , (6) 
where the y H : y c is 1 :0.25 are the gyromagnetic ratio of the nuclei. We describe the various stages 
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FIG. 2: The pulse sequence for the implementation of PPS and T operator. The pulses are applied from 
left to right. All the black narrow pulses have a flip angle of 90° and the white broad pulse have a flip angle 
of 180°. The flip angle of the other pulses are written below them. The phases of the pulses are written 
above them.Gj represent z-gradients which are used to kill unwanted transverse magnetization. The time 
period of jj implies that the system is evolved only under the scalar coupling Hamiltonian and the evolution 
under the chemical shifts are refocused by the n pulses. The flip angle = | , | or | depending upon a and 
®=2fs. 

of the experimental implementation in the following paragraphs. 

For a two-qubit system, there are 4 basis states: |00), |01), |10) and |1 1). We choose the target 
state to be |11). If \s) is an equal superposition of the basis states then a = (s\ll) = 1/2. But 
to simulate the Grover's algorithm for large database, we want a to be small. That we achieve 
by letting \s) to be an unequal superposition. We first create the |00) PPS by the use of spatial 
averaging lIlOn .rFig. |2j. A |00) PPS has a unit population in the |00) state and zero population in 
|01), 1 10) and |11) states. Then we apply a V pulse on it. We have 



100) 















cos z I — I |00) + sin ( - 1 cos | - 1 101> 
+ s,n(|jcos(|)|10) + sin 2 (|||l!) 



(7) 



Thus a = (s\U) = sin 2 (f ). By choosing = f , f and §, we achieve a = 0.146, 0.067, 0.030 
respectively. Just to compare, if \s) is an equal superposition of /V basis states then a = 1 / and 
these values of a correspond to N x 47, 223, 1111 respectively so that we need n = log 2 N « 
6, 8, 10 qubits respectively to represent all basis states. However, by choosing \s) to be an unequal 
superposition, a two-qubit system becomes sufficient to simulate large databases. 

The next step in the implementation of the algorithms is the application of the Q 2 IT operator. 
In our case, we assume that there are no errors in I t trasnformation, i.e. <p = n. Since |s) = y |OO) 



and \t) = |11), we have 

@ 2 = ® y I^ y In®yf m ®lhi , (8) 
T = ® y I^®;in®y4o®llu ■ (9) 

Note that in case of no errors in Too transformation, we have <f> = n. Fig. |2] contains the pulse 
programme for the implementation of the T operator. The EL and In operators are selective phase 
rotations of |00) and |1 1) states respectively i.e. 

f 

|00) + |01) + |10) + |11> A e'*|00> + |01> + |10> + |ll>, (10) 
|00) + |01) + |10> + |11) |00) + |01) + |10)-|11). (11) 

Therefore in NMR, the 7^ and In operators are implemented by evolution under 

C = exp(/ z 1 + / z 2 + 2/^) (12) 
Iu = <xp(-ll-% + 2ltf) (13) 

respectively. Following [5], for 1^, the evolution under l\ and I 2 are implemented by composite 
z-rotation pulses like [| [j] _ . The evolution under the 1l\l 2 z is implemented by evolving the 
system under the scalar coupling Hamiltonian only, for a time period of 4>l2nJ . The I q operator is 
applied by (a) reversing the order of application of pulses and evolution, (b) flipping the phase of 
the centre pulse of the composite z rotation by n and (c) changing the evolution time from cp/lnJ to 
(An-<p)/2nJ for /if. The application of In involves evolution of the system under the Hamiltonian 
given by Eq. [131 This is similar to the Hamiltonian evolution for I^ , the only difference being the 
negative sign before 1} and I 2 . This implies that the phase of the centre pulse in the composite z 
rotation is changed by n. Moreover, as errors are not introduced in the In operator (i.e. (p = n), 
the central pulse in the composite z rotation has a flip angle of | and the time of evolution is 1/2J 
in both the cases of implementation. 

After the implementation of the algorithm, the final state is measured. In this case only the 
diagonal elements of the final density matrix (population spectra) is required to be measured. This 
is done by collecting the data after applying a gradient to kill the off-diagonal elements followed 
by a 90° pulse. The diagonal elements of the final density matrix are reconstructed from the 
population spectrum. 

Fourteen iterations of the original and the modified algorithms were implemented for three 
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FIG. 3: The plot of 'probability of the marked state' calculated from the measurement of the population 
(see text) of the final state for each iteration of the 'original' (depicted by open circles) and the 'modified' 
(depicted by crosses) search algorithm. The line joining the points are a guide to the eye. As the value of 
decreases a also decreases i.e the size of database increases, (a) = | and no error in I s or I t . Both the 
'original' and the 'modified' search algorithm yield the marked state with high probability, (b) = 5 and 
a 10% error in I s . The 'original' algorithm cannot amplify the amplitude of the marked state to the desired 
level while the modified algorithm is able to search the marked state with a much higher probebility. (c) 
= | and 10% error in I s . The 'original' algorithm has failed completely in searching the marked state, 
while the 'modified' algorithm succeeds. 



different values of initial probability amplitude of the marked state i.e. = |, | and |. For the 
case = j, no error has been introduced in 7 o operator (i.e. (f> = n) which implies that the 
phase matching condition is satisfied in this case. It can be seen that both 'original' and 'modified' 
algorithm behaves almost similarly i.e. they find the marked state with almost the same periodicity 
(Fig. |3(a)[ ). In the next case (Fig. |3(b)[ ), the value of as chosen to be | to make a smaller and 
an error of 10% was introduced in Too (i.e. 4> = 0.9n) so that the phase matching condition is 
violated. We see that in this case, the original search algorithms starts to fail while the 'modified' 
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FIG. 4: Bar-plot of the diagonal elements of the output density matrix, (a) after PPS (c) After application 
of six steps of 'modified' search algorithm for = | and error in l is 10%. The marked state |11) has a 
high probability signifying the sucess of the algorithm. 

algorithm obtains the searched state with a high probability (~ 80 %). The original algorithm 
cannot amplify the amplitude of the marked state so as to definitely distinguish it and therefore 
the solution is not reached. Finally, the algorithms were implemented for = | and 10% error in 
Too operator (Fig. |3(c)[ ). In this case a is very small (simulating a system of about 10 qubits), and 



therefore the 'phase matching' condition is violated even more strongly. It can be seen that in this 
case, the 'original' algorithms totally fails in reaching the solution but the 'modified' algorithm 
succeeds. For completeness, the diagonal elements of the tomographed density matrix for the case 
of Fig. |3(c)| i.e. = | and 10% error in 7 o are plotted in Fig 4. This confirms the success of the 
'modified' algorithm of Tulsi [|70. 



In conclusion, we have implemented the 'modified' quantum search algorithm by Tulsi [|7|] 
and have experimentally validated his claim that his algorithm is robust to errors in U G operator 
as compared to the original search algorithm. We have shown that small errors can be fatal for 
searching larger databases using Grover's algorithm while the 'modified' search algorithm is ro- 
bust. We have experimentally simulated the behaviour of the algorithms in large database on a 
2-qubit NMR quantum information processor. Quantum computers when fully operational will be 
dealing with real world problems requiring large systems. This experiment, besides providing a 
validation for an important theoretical prediction, will help in providing impetus to future work on 
the study of existing algorithms in large real world systems. 
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